Scaling and Eigenmode Tests of the Improved Fat Clover Action 

Mark Stephenson and Carleton DeTar 
Department of Physics, University of Utah, Salt Lake City, UT 84112, USA 

Thomas DeGrand and Anna Hasenfratz 
Department of Physics, University of Colorado, Boulder, CO 80309, USA 

(February 7, 2008) 

We test a recently proposed improved lattice-fermion action, the fat Unk clover action, exam- 
ining indicators of pathological small-quark-mass lattice artifacts ( "exceptional configurations" ) on 
quenched lattices of spacing 0.12 fm and studying scaling properties of the light hadron spectrum 
for lattice spacing a = 0.09 and 0.16 fm. We show that the action apparently has fewer problems 
with pathological lattice artifacts than the conventional nonperturbatively improved clover action 
f^ ' and its spectrum scales just as well. 
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I. INTRODUCTION 



The goal of lattice ferniion improvement schemes is to increase the effectiveness of computer algorithms by signif- 
icantly reducing lattice artifacts at economically coarse lattice spacings, particularly in the range a < 0.15 fm. The 
Cf^ ' clover action with the clover coefficient chosen nonperturbatively [|] (NPCA) has been found to give a substantially 
K*" ', better scaling of the hadron spectrum over the range [0.06,0.16] fm than the conventional Wilson fermion action 
^^ (WA) Ig,^. However, as the lattice spacing is increased over this range, the NPCA becomes increasingly sensitive 

^i^ \ to local fluctuations in the gauge configurations that produce unwanted artifact singularities at small, positive quark 
mass, the so-called "exceptional configurations." To avoid them, one must keep the quark mass artificially high. At 
increasingly coarse lattice spacing, the lower bound on "safe" quark masses rises, making an extrapolation to physical 
^\ , quark masses increasingly problematical. Difficulties with such "exceptional" configurations are overcome in a variety 
0^ of approaches, including pole shifting W and schemes which implement an exact chiral symmetry on the lattice, the 
^^ : overlap formalism 0, and the domain wall approach 0. Such methods are computationally expensive. Here, we 
^ consider a new (approximate) improvement scheme, the "fat clover" action. 

I ' The "fat clover" action (FCA), proposed by DeGrand, Hasenfratz, and Kovacs m, couples the standard clover 

^^' action to a locally smoothed gauge field. Smoothing of the gauge fields is achieved through a series of APE blocking 
f-H ' steps S. It is intended that the number of blocking steps remains fixed as the continuum limit is approached. 
• • , Thus the fermion-gauge coupling is modified at a scale that is a fixed multiple of the cutoff, and the correct local 
'^ ' action is recovered in the limit. Smoothing has a number of beneficial effects: lattice artifacts are suppressed, chiral 
K^ , properties are improved, and the renormalization of a variety of lattice quantities, such as the local-vector-current 
5^ ■ and axial- vector-current renormalization constants Zy and Za, is small mM- 

Cu I Here we examine two variants of fat link actions: one using the tree level value for the strength of the clover term 

(TFCA) and one using an optimized value (OFCA). (See Table | for a guide to our abbreviations.) 

We first study the distribution of real eigenmodes of the lattice Dirac operator in an ensemble of gauge configurations, 
with particular attention to the "exceptional" eigenvalues at positive quark mass. We show that the FCA has improved 
chiral properties, in the sense that the spread of leading near zero modes is narrower than with the NPCA at a lattice 
spacing of 0.12 fm. Thus the "safe" lower bound on quark masses is lower with the FCA. We then present a study of 
fluctuations in the pion correlator as a qualitative test of the suppression of exceptional configurations in NPCA and 
TFCA. 

Finally we perform a small scaling test of the quenched hadron spectrum for the TFCA and OFCA. Scaling is tested 
as the lattice constant is varied from 0.092 fm to 0.164 fm, choosing some fixed value of the quark mass, such that 
either rn^ = 2.5(T, where a is the string tension, or m-^/rrip « 0.7 H]. Included for comparison are the corresponding 
results for the standard Wilson action (WA), nonperturbative clover action (NPCA), and the standard Wilson action 
on a fat link background (FWA). In this scaling study we pay careful attention to the elimination of a variety of 
sources of systematic error: we scale all quantities with physical dimensions and fix the lattice dimensions in physical 

units. 

In Sec. II we give details of the FCA and our computational method. In Sec. HI we present results of our lattice 



simulations. We conclude with Sec. [V 



II. COMPUTATIONAL METHOD 
A. Fat Clover Action 

The fat clover action (FCA) is the usual clover action 

Ssw = ipM{K)ip ^ Sw - Kcsw ^ ip{x)af_,^Pf_,^ip{x) (1) 

where Sw is the Wilson action, and Pp,y is the standard "clover" expression for the field strength tensor F^,y, except 
that all link variables in Sw a-nd Pfj_i, are fat links. The fat link is constructed from the original ("thin") links with a 
series of APE blocking steps |^ . A single step creates a new gauge configuration with each gauge link replaced by a 
weighted sum of the link and its staples, followed by a projection back to SU(3). Explicitly, each link Ux f_i is replaced 

by 



T4,^ = V 



(1 - c)Ux,f^ ~^ qZ^ [Ux,i^Ux+i,,^i,U^^^^ + U^_^i^Ux-u,t^Ux-u+p,,i' 



y^fi. 



(2) 



where V represents the projection, which in SU(3) chooses the unique group element U of SU(3) maximizing Tv{UV^). 
(Replacement occurs after all smooth links for the lattice are computed.) This process is repeated TV = 10 times in 
our study with a coefficient c = 0.45. Such a choice was found to give good stability in instanton size and placement 
during smoothing in both SU(2) and SU(3) [|l^,0. We have not explored other choices extensively, but suspect 
considerable latitude is permitted in the choice of c and N . After smoothing, the mean plaquette of the fat link is 
close to 1: For 6/g^ = 5.7, the mean plaquette is Tr[/p;aq/3 — 0.985 and for 6/g^ — 6.0 it is 0.994. 

Fat link clover actions, like all clover actions, have no 0{a) nor 0{g^a) corrections as long as the clover coefficient 
csw is taken to approach unity in the 5^0 limit. In perturbation theory fl^ the effect of the fat link is to multiply 
the usual thin link quark-gluon vertices by a form factor (1 — ^a^q^)^ with q^ = (2/a) sin(g^a/2), where q is the gluon 
momentum. This easily accounts for the observed considerable reduction in additive quark mass renormalization and 
finite renormalization constants for vector and axial currents close to unity. 

The clover coefficient csw is a priori unspecified. We expect, based on the near unit value of the mean plaquette, 
that an optimal value for simulations would be close to the tree-level value of unity, and that is what one would choose 
in a tadpole-based improvement program [ |l3| . We chose the value of csw using the approach of DeGrand, Hasenfratz 
and Kovacs, based on the position of real eigenmodes of the Dirac operator [Q. 

To locate the real eigenmodes, for each configuration we calculated a noisy estimator of the expectation value 

(V^TsV') 000 

A(«)=^75Af -1(^)7? (3) 

at a closely spaced series of real values /t, where 77 is an arbitrary random vector, held constant for the scan over k. A 
real eigenvalue of Af (k) appears as a pole in A{k), provided the corresponding eigenvector has nonzero overlap with 
the vector 77. The quantity fjM~^{K)r] would also diverge at a real eigenvalue of M{k), but since the real eigenmodes 
of M{k) are also close to being eigenmodes of 75 with eigenvalue ±1, the factor 75 helps to distinguish them. 

In the continuum limit the only real eigenvalues of M(k) are chiral zero modes, occurring at Kc = 1/8, i.e. zero bare 
quark mass. If the lattice Dirac operator is not chiral, the real modes are spread around Kc (defined as the value of k 
where m'^{K) extrapolates to zero), and can also be shifted from the chiral Kc = 1/8 value. Such real eigenmodes are 
undesirable lattice artifacts that prevent lattice simulations at small quark mass. ||I^. Configurations that produce 
them are called "exceptional" , although the problem really lies with the choice of fermion action and not with the 
gauge configuration itself. Actions with improved chiral properties have real eigenmodes that cluster more closely 
around Kc. For such actions it should be possible to study lower quark masses without encountering difficulties with 
exceptional configurations, or to carry out a simulation on a coarser lattice at the same quark mass. 

For a sufficiently high degree of fattening one can optimize the clover coefficient by minimizing the spread of the 
real eigenmodes on Monte Carlo-generated configurations. In this work, however, we choose a simpler approach. We 
generate a series of artificial lattice instantons of varying size a < tq < 3a on 8 lattices and of size 2a < rg < 6a on 
12^^ lattices. Instanton studies at 6/(7^ = 5.7 and 6/(7^ = 6.0 predicted instanton sizes in this range [Q. For each such 
gau ge c onfiguration, we examined the position of the resulting fermion real eigenmode. Results are shown in Tables O 
and [I] and Fig. |l|. For large tq the near-zero modes are quite close to zero quark mass. As rp drops below the lattice 



cutoff, the would-be zero mode moves toward negative quark mass. The trajectory of real eigenmodes is altered by 



adjusting the clover coefRcient csw- We find that with csw = 1-2 the variation in pole position is minimized for 
instanton sizes in the range rg > a. With csw = 1-1 the variation in pole position is minimized for rg > 2a. Since our 
scaling test considers lattices over a range of spacings varying by a factor of 2, we choose csw = 1-2 for our coarsest 
lattice spacing and csw = 1-1 for our finest lattice spacing to assure scaling consistency. We call the action with 
tuned clover coefficient the "optimized fat clover action" (OFCA). To test sensitivity to this choice, we also present 
results with the tree-level choice csw = 1 (TFCA). 

B. Computational Parameters 

There are two parts to this study: an analysis of the distribution of low energy real eigenvalues and an analysis 
of the scaling of the spectrum. For studies of the distribution of near-zero eigenvalues at 6/5^ = 5.85 we analyzed 
100 configurations of size 10''', and for the companion spectrum study to relate the pion mass to the k, value, 20 
configurations of size 12^ x 48. For the spectrum scaling studies we have worked with two ensembles of quenched 
gauge configurations generated with the conventional one-plaquette Wilson action: 120 configurations of size 8^ x 24 
at 6/g^ = 5.7 and 100 of size 16'^ x 48 at 6.0, corresponding to lattice spacing a = 0.164 and 0.092 fm, respectively, 
based on recent measurements of the string tension for this action in lattice units [Q and the choice \/a = 468 MeV 

The lattice dimensions for the scaling study were chosen to keep an approximately constant physical volume, so as 
to avoid inconsistent finite size effects. To allow tuning of the quark mass, either by fixing the pion mass in terms of 
the string tension or in terms of the rho mass, we calculate the spectrum for (typically) three neighboring k values 
selected so that the desired dimensionless ratios can be reached by interpolation. Quark propagators were generated 
from a fixed Gaussian (shell- model) source with standard deviation 2a for the coarsest lattice and 4a for the finest. 
Our scaling tests were computed at three values of the clover coefficient csw for each gauge coupling: csw = 1-1 or 
1.2, as noted above, for the OFCA, csw = 1 for the tree-level fat clover action (TFCA), and csw = for the fat 
Wilson action (FWA) to test the relative merits of smoothing and reducing the 0(a) errors in the action. 

III. RESULTS 
A. Fermion eigenmodes 

For each of the actions in our study we determined the distribution of the real eigenmodes on a set of 100 10^ gauge 
configurations, generated with the conventional single plaquette action at 6/g^ = 5.85. We determined the probability 
distribution P{m'^) of the leading pole (i. e. the eigenvalue corresponding to the largest quark mass) for the various 
actions in our study. 

Note that this statistic is different from the eigenvalue histograms of Ref . H , where all the low energy eigenmodes 
were included to study the spread of the physical modes. Since here we consider only the pole corresponding to the 
largest mass on each configuration, these plots are indicative of the exceptional configurations. 

To compare leading-pole distributions from the various actions, we converted the k values to m^ values, by measuring 
the hadron spectrum with the same action on a set of 20 12'^ x 48 quenched configurations. Results are shown in 
Fig. g. Pion masses used for constructing the linear scale conversion, rn^ = a/n + b are given in Tables [V and|^. The 



bin widths in this figure are variable, since the poles were located by scanning at the same constant increment in k for 
all actions. We took this interval as the resolution of the pole location. The corresponding interval in m^, however, 
varied from action to action. The bin heights are scaled so that the probability distributions all have unit total area. 
The TFCA and OFCA actions clearly produce distributions that are more sharply clustered around m^ = 0. The 
peak for the OFCA appears at a nonzero bin in m^ , but we estimate a combined systematic and statistical error (one 
sigma) of one bin width arising from the conversion from k to rn^ near zero pion mass. 

To put these results in another perspective, we have also measured the ratio rriT, /mp for these configurations. For 
our sample of gauge configurations aX (3 — 5.85 the NPCA encounters its first pole at ra.,^/mp = 0.56(4), the TFCA, 
at about 0.37(6) and the OFCA at about 0.42(5). 

B. Pion Correlators 

While a pole in the quark propagator at k < Kc is the most precise indicator of an exceptional configuration, 
fluctuations in a hadron correlator, e.g. the pion correlator, give a qualitative indication. Poles in the quark propagator 



typically have residues with concentrated support in Euclidean space-time. For instanton-induced poles, localization 
of the residue comes from localization of the zero eigenmode. Since the pion propagator is the gauge-invariant square 
of the quark propagator, a nearby exceptional pole typically contributes a strong localized fluctuation in the pion 
correlator— in some cases inducing a "W" shape in a semilog plot. 

Starting from a common sample of 80 quenched single-plaquette 8^ x 24 gauge configurations at 6/5^ — 5.7, we 
compute the pion correlator for the NPCA and TFCA for a range of quark masses. We then consider two measures 
of fluctuations: (1) the noise to signal ratio of the correlator at a fixed time and (2) the number of outliers, based on 
a correlated chi square measure. 

The simplest measure of fluctuations is the noise to signal ratio in a correlator, namely the ratio of the standard 
deviation of a correlator a{t) to the value of the correlator, c(i), as a function of distance and/or quark (or pion) 
mass. The pseudoscalar correlator is the most useful one to look at, since simple theoretical arguments suggest that 
a{t)/c{t) should be roughly independent of the quark mass. In Fig. ^ this ratio at i = 10, namely tT(10)/c(10), is 
plotted over a range of mT^/mp for the NPCA and TFCA. It is clear that for 0.5 < mT^/mp < 0.7, fluctuations are 
dramatically reduced with the TFCA. 

We denote the correlator on the ith configuration by Ci (t) and its mean over the sample of configurations by c(t) . 
For the outlier test we start by constructing the usual covariance matrix Wt.f, based on the observed fluctuations, and 
its inverse Wt.f. The correlated chi square measure for configuration i is then 

xi-J2^cm~m]hit')-cit')]w,,,, (4) 

for Nt degrees of freedom. The corresponding confidence level is then used to determine the strength of deviation 
from the mean. For the sample of 80 configurations we treated the two time intervals [0, Nt/2] and [Nt/2 + l,Nt — 1] 
separately. A configuration was deemed exceptional, if the confidence level determined on either interval was less than 
10~^, a somewhat arbitrary value that was chosen to correspond to strongly discernible deviations, many of them with 
the "W" shape, characteristic of an exceptional configuration. Note that a Gaussian normal fiuctuation at this level 
in a sample of 80 would be expected only once in about 10^ trials. The process of identifying outliers was carried out 
iteratively, in each pass removing the outliers from the sample as they were identified, until none remained. Results 



are summarized in Tables VI and VII. Based on this measure, if we were to insist on no more than one exceptional 
configuration in a sample of this size at /3 = 5.7, the NPCA would be restricted approximately to rriTr/mp > 0.8 and 
the TFCA to m^/nip > 0.6. 

C. Spectrum 

Correlators for the zero momentum pion, rho, nucleon, and A were fit to single exponential forms, minimizing the 
correlated x^. Care is needed to prevent biases arising from the choice of the fitting range, particularly from the 
choice of minimum time imin. We used two methods to test for bias: In fitting correlators for a set of closely spaced k 
values we selected (1) the smallest tmin giving a minimum CL for all k's greater than 0.05 and an average CL greater 
than 0.1; and (2) the t,nin for which the product of CL and the number of degrees of freedom (df) is maximum, a 
rather ad hoc rule of thumb [2^. As a rule both methods gave the same tmm- Where a different value was obtained, 
we determined that the variation in mass value was within the statistical errors of the fits. 

We found that tmin for the 8^^ x 24 lattice was generally half the value on the 16^ x 48 lattice. In the few cases 
in which it was not (the rho meson for OFCA and TFCA), we verified that, had we enforced this further condition, 
the central mass value would have shifted by less than 1%, an amount smaller than the error in the observed scaling 
violation. Our fitting range is then approximately constant in physical units, and our results are therefore free of bias 
from this source. 

We also checked the single exponential fits against two-exp onen tial fits and verified that the results were stable 



within statistical errors. Results of the fits are shown in Table VIII. 

We consider two alternatives for fixing the quark mass: (1) fixing m^/a — 2.5 and (2) fixing 111^^ /nip — 0.7. These 
values were chosen to correspond to each other, approximately. Since variations in the strength of the clover term 
changes the A^ — A and n — p mass splittings, so can change the pi to rho mass ratio at fixed physical quark mass, the 
former method is preferable. We present the second, more popular, method to allow comparison with other work. 

Table [X and Figs. |j and |g show the masses of particles and their ratios at mj = 2.5cr and the scaling violations from 
6/g^ = 5.7 to 6/g^ = 6.0. For comparison, conventional Wilson data (WA) from recent calculations are shown. The 
WA 6/g^ — 5.7 values are interpolated from raw data given in Rcf. |gl[. The WA 6/g^ = 6.0 values are interpolated 



from raw data given in Refs. p3,E3]. The value rn^ = 2.5cr is slightly outside the range of values given in Ref. |22|, so 



to avoid extrapolation, those data were supplemented by data given in Ref. ||23||; however, extrapolation from the data 



of Ref. 12^ alone gives the same results. The 6/g^ = 5.7 NPCA values were interpolated from unpublished values 
provided by Heller [|||. The 6/5^ = 6.0 NPCA values are from Ref. [||. 

The mass of the rho meson is seen to be a sensitive indicator of scaling violations. The scaling violation in the rho 
mass is reduced from approximately 7% for the WA to less than 2% for the OFCA and TFCA. Fattening the WA 
does not improve scaling with any significance. Thus smoothing of the gauge fields alone does not improve scaling. 
This is because the Wilson action, with either thin or fat links, has 0(a) lattice artifacts which are removed by the 
addition of the clover term. (A similar behavior for a fat link action without a clover term, in that case, a hypercubic 
action, was also seen in Ref. p5|.) 

Scaling violations of the nucleon mass are statistically consistent for all the actions considered. Scaling violations 
of the delta mass are improved from roughly 8% for the FWA to less than 4% for the TFCA and OFCA. This also 
is found for the alternative fits of the delta masses, suggesting that the scaling improvement seen is real, despite 
potentially large systematic errors in the fits of the delta masses. 

We extrapolated the mass values to zero lattice spacing, forcing a common extrapolated mass for all actions. Our 
extrapolation is linear in a for the WA and FWA and hnear in a^ for the OFCA, TFCA and NPCA. Results are 
plotted in Figs. ^ and Q It is clear that none of the actions completely remove scaling violations in the nucleon or 
delta mass, but that all of the clover actions show smaller violations than the Wilson actions. 

For the second approach we adjust quark masses so as to fix the ratio niTj/nip = 0.7. Table ^ shows the masses 
of particles and their ratios, and the scaling violations from 6/g^ = 5.7 to 6/g^ = 6.0. For comparison, conventional 
Wilson data from recent calculations are shown. The Wilson 6/g^ — 5.7 values are interpolated from raw data given in 
Ref. pw. The Wilson 6/g^ = 6.0 values are interpolated from raw data given in Ref. l22|. Also shown are data for the 
TFCA and NPCA. For the NPCA the 6/g^ = 5.7 values are given in Ref. [D, already interpolated to mT,/mp = 0.7. 
The values at 6/3^ = 6.0 are from Ref. |^. 

The mass of the rho again is seen to be a sensitive indicator of scaling violations. For the pair of lattice spacings 
used, the scaling violation of the rho mass is reduced from approximately 13% for the conventional WA to less than 
2% for the OFCA or the TFCA. There is no significant scaling improvement of the rho mass for the FWA compared 
with the WA. The scaling violations of the rho are compounded with those required by the fixing of the pion to rho 
mass ratio to the same constant for all cases. The quark masses are forced to values at &/ g^ — 5.7 and 6/g^ = 6.0 
such that the scaling violation of the pion mass equals that of the rho mass. This roughly doubles the total rho 
scaling violation at fixed m-^/nip, compared with that for fixed m-^, for cases where scaling violations in the pi-rho 
mass splittings are large, as will be shown to be the case for the Wilson actions. 

Scaling violations of the nucleon mass also are compounded with those required by the fixing of the pion to rho 
mass ratio to the same constant for all cases. The value of the nucleon mass is lowered similarly to the pion mass by 
the clover (magnetic moment) term, so forcing the scaling violation of the pion mass to equal that of the rho mass also 
forces the nucleon mass to acquire a similar scaling violation. At m^/mp = 0.7, scaling violations of the nucleon mass 
are reduced from approximately 10% for the WA and FWA to approximately 5% for the TFCA to approximately 2% 
or lower for the OFCA. 

Scaling violations of the delta mass are improved from roughly 14% for the FWA to less than 5% for the OFCA. 



Mass splittings are given in Table XI at m^ = 2.5cr. Phenomenological values of the differences of the squared masses 
of vector and pseudoscalar particles are almost equal for different quark flavors, so the unrealistically high quark mass 
used in the lattice calculations should not matter much. Similarly, the difference of the masses of spin -^ and spin 
I baryons are comparable for different flavors. Scaling violations of the mass splittings are reduced substantially for 
the clover actions compared with the Wilson actions. Also there is rough agreement with experimental values for the 
clover actions. The experimental values are {m? — mj) — 5.7 x 10^ MeV^ and toa ~ "rriN — 294 MeV. Using ^Ja — 468 
MeV, we have (m^ — m1)/a = 2.6 and {ttin — mA)/\/a- = 0.63. 

In summary, the proposed fat clover actions with either tree-level clover coefficient (TFCA) or optimized clover coef- 
ficient (OFCA), have greatly improved scaling properties compared with the Wilson action. The scaling improvement 
with either of these actions is comparable to that of the NPCA. Scaling tests of the FWA show that smoothing by 
itself does not improve scaling. The combination of the clover term with a fat link allows one to reach smaller values 
of the pseudoscalar mass than is possible with the NPCA, with an apparently equivalent level of scaling violations. 

IV. CONCLUSIONS 

Using a common set of gauge configurations, we have carried out a systematic study of the distribution of leading 
near-zero eigenvalues and the scaling of the light hadron spectrum and for a variety of fermion actions on quenched 
lattices with lattice spacing in the range [0.09,0.16] fm. Actions included in this study are the conventional Wilson 
action (WA), Wilson action on fat gauge links (FWA), clover action with a non-perturbatively tuned clover coefficient 



(NPCA), tree-level "fat" clover action (TFCA) and optimized fat clover action (OFCA). All of the clover actions 
show better scaling behavior than any of the Wilson actions-an entirely expected result, since adding the clover term 
with the correct coefficient converts an 0{a) action into an 0{a^) one. Based on an analysis of pole positions, we 
have found that the fat clover actions OFCA and TFCA exhibit chiral properties superior to the NPCA, WA, and 
FWA. A further analysis of fluctuations in the pion correlator shows that the TFCA is far less noisy than the NPCA, 
a further circumstantial indication of a suppression of exceptional configurations. 
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TABLE L Guide to abbreviations for the fermion actions in this study. 

FGA Clover action on fat gauge links (either optimized or tree- level). 

OFCA Clover action on fat gauge links with optimized clover coefficient. 

TFCA Clover action on fat gauge links with tree-level clover coefficient. 

NPCA Conventional clover action with non-perturbatively tuned clover coefficient. 

FWA Wilson action on fat gauge links. 

WA Conventional Wilson action (r = 1). 
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for small instantons 
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1.2 
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1.8 


1.00 
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-0.003 
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1.25 


-0.083 
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0.012 
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1.50 
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0.021 
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0.146 
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-0.030 
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0.025 
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0.194 


2.25 


-0.012 


0.006 


0.025 




0.045 


0.088 


0.159 


2.50 


-0.009 


0.007 


0.023 




0.041 


0.078 


0.138 


2.75 


-0.007 


0.007 


0.021 




0.037 


0.070 


0.122 


3.00 


-0.006 


0.006 


0.020 




0.034 


0.064 


0.111 



TABLE III. Pole location rUpa for larger instantons 



p/a 



1.0 



csw 

1.1 



1.2 



2.00 
2.50 
3.00 
3.50 
4.00 
4.50 
5.00 
5.50 
6.00 



-0.019 
-0.008 
-0.004 
-0.002 
-0.002 
-0.002 
-0.002 
-0.002 
-0.003 



0.003 
0.005 
0.005 
0.005 
0.004 
0.003 
0.003 
0.002 
0.002 



0.024 
0.020 
0.015 
0.012 
0.010 
0.009 
0.007 
0.007 
0.006 



TABLE IV. Pion masses (smeared-local channel) in lattice units vs k for a variety of actions on quenched 12'' x 48 lattices 
at 6/g2 = 5.85. 





OFCA 






TFCA 






NPCA 




0.1200 




0.621(6) 


0.1220 




0.539(4) 


0.1300 




0.733(7) 


0.1210 




0.547(6) 


0.1230 




0.463(4) 


0.1310 




0.645(7) 


0.1220 




0.465(7) 


0.1240 




0.378(5) 


0.1320 




0.547(8) 


0.1230 




0.370(11) 


0.1245 




0.319(5) 


0.1330 




0.432(9) 


0.1235 




0.316(17) 


0.1250 




0.260(12) 


0.1335 




0.364(10) 



TABLE V. Continuation of table | iv| 





FWA 






WA 




0.1265 




0.515(6) 


0.1560 




0.511(9) 


0.1275 




0.455(6) 


0.1570 




0.457(10) 


0.1285 




0.389(7) 


0.1580 




0.400(12) 


0.1295 




0.309(10) 


0.1590 




0.336(16) 


0.1300 




0.261(13) 


0.1595 




0.298(18) 



TABLE VI. Outliers at CL < 10"'^ for the NPCA in a sample of 80. 



niTr/mp 



N 



0.125 

0.127 

0.129 

0.1295 

0.130 

0.1302 

0.1303 



0.84(1) 
0.78(1) 
0.68(1) 
0.64(1) 
0.48(2) 
0.50(2) 
0.51(3) 




4 
5 
7 
12 
15 
14 



TABLE VIL Outliers at CL < 10"^ for the TFCA in a sample of 80. 



m-^/mp 



N 



0.121 
0.122 
0.123 
0.124 
0.125 
0.126 



0.80(1) 
0.75(1) 
0.70(2) 
0.62(2) 
0.50(2) 
0.42(4) 



TABLE VIII. Summary of hadron masses and bootstrap errors for various fat-link actions: OFCA {csw = 1.1, 1-2), TFCA 
(csw = 1.0), and FWA (csw = 0.0). 



csw 


6// 


K 


mTrffl 


nipa 


TUN a 


mAfl 


1.1 
1.1 
1.1 


6.0 
6.0 
6.0 


0.1225 
0.1230 
0.1235 


0.3985(21) 
0.3549(23) 
0.3072(25) 


0.539(5) 
0.511(6) 
0.483(7) 


0.807(7) 
0.752(8) 
0.691(9) 


0.899(11) 
0.865(14) 
0.832(17) 


1.2 
1.2 


5.7 

5.7 


0.1200 
0.1220 


0.731(4) 
0.569(5) 


0.968(7) 
0.882(10) 


1.411(17) 
1.223(24) 


1.58(2) 
1.44(3) 


1.0 
1.0 
1.0 


6.0 
6.0 
6.0 


0.1225 
0.1230 
0.1235 


0.4293(20) 
0.3891(22) 
0.3460(25) 


0.553(5) 
0.526(5) 
0.498(6) 


0.841(7) 
0.791(7) 
0.737(8) 


0.920(10) 
0.882(12) 
0.848(14) 


1.0 
1.0 

1.0 


5.7 
5.7 
5.7 


0.1200 
0.1220 
0.1245 


0.816(4) 
0.682(4) 
0.478(6) 


0.992(5) 
0.907(6) 
0.802(10) 


1.496(15) 
1.335(16) 
1.095(28) 


1.630(19) 
1.503(23) 
1.341(29) 


0.0 
0.0 
0.0 


6.0 
6.0 
6.0 


0.1270 
0.1280 
0.1290 


0.355(3) 
0.280(4) 
0.186(5) 


0.455(5) 
0.410(5) 
0.369(9) 


0.725(6) 
0.631(7) 
0.519(13) 


0.802(13) 
0.743(13) 
0.678(22) 


0.0 
0.0 
0.0 


5.7 
5.7 
5.7 


0.1280 
0.1310 
0.1330 


0.666(5) 
0.501(7) 
0.435(9) 


0.780(6) 
0.681(8) 
0.650(8) 


1.263(18) 
1.073(22) 
0.993(22) 


1.359(17) 
1.204(23) 
1.151(27) 



TABLE IX. Scaling summary at m^ = 2.5ct 



6/ff^ 


mT,/mp 


ruN Imp 


mp|^/a 


mjv/^ 


•yn/^l^/a 








WA 






5.7 


0.815(5) 


1.596(11) 


1.929(21) 


3.075(36) 


- 


6.0 


0.763(2) 


1.538(6) 


2.072(10) 


3.186(17) 


- 




-6.4(0.6)% 


-3.6(0.8)% 


+7.4(1.3)% 
FWA 


+3.6(1.3)% 




5.7 


0.823(11) 


1.604(30) 


1.921(26) 


3.082(59) 


3.362(60) 


6.0 


0.768(11) 


1.587(21) 


2.055(23) 


3.261(31) 


3.631(62) 




-6.6(1.7)% 


-1.1(2.3)% 


+7.0(1.9)% 
TFCA 


+5.8(2.3)% 


+8.0(2.7)% 


5.7 


0.704(8) 


1.439(24) 


2.246(29) 


3.233(61) 


3.733(75) 


6.0 


0.695(10) 


1.480(23) 


2.275(29) 


3.366(38) 


3.872(65) 




-1.3(1.8)% 


+2.8(2.4)% 


+1.3(1.8)% 
OFCA 


+4.1(2.3)% 


+3.7(2.7)% 


5.7 


0.678(9) 


1.408(29) 


2.334(34) 


3.285(66) 


3.814(80) 


6.0 


0.684(10) 


1.462(24) 


2.310(30) 


3.384(39) 


3.922(68) 




+0.9(2.0)% 


+3.8(2.7)% 


-1.0(1.9)% 
NPCA 


+3.0(2.4)% 


+2.8(2.8)% 


5.7 


0.671(8) 


1.432(20) 


2.356(30) 


3.375(52) 


3.854(87) 


6.0 


0.674(2) 


1.488(9) 


2.345(12) 


3.490(23) 


- 




+0.4(1.2)% 


+3.9(1.6)% 


-0.5(1.4)% 


+3.4(1.7)% 


- 



6/5^ 



rati I "trip 



TABLE X. Scaling summary at rxiT^ jnip — 0.7 



mNJ^fo 



mA/\f5 



5.7 
6.0 



5.7 
6.0 



5.7 
6.0 



5.7 
6.0 



5.7 
6.0 





WA 






1.554(16) 


1.698(19) 


2.638(35) 


— 


1.514(8) 


1.924(10) 


2.913(17) 


- 


-2.6(1.1)% 


+13.3(1.4)% 
FWA 


+10.4(1.6)% 




1.547(38) 


1.710(27) 


2.646(62) 


3.024(71) 


1.548(24) 


1.905(24) 


2.952(33) 


3.463(61) 


+0.1(2.9)% 


+11.4(2.3)% 
TFCA 


+11.6(2.9)% 


+14.5(3.4)% 


1.436(24) 


2.247(29) 


3.235(61) 


3.735(75) 


1.475(22) 


2.289(28) 


3.393(37) 


3.889(64) 


+2.7(2.3)% 


+1.9(1.8)% 

OFCA 


+4.9(2.3)% 


+4.1(2.7)% 


1.422(26) 


2.383(33) 


3.393(64) 


3.893(76) 


1.475(22) 


2.349(29) 


3.465(38) 


3.968(65) 


+3.7(2.4)% 


-1.4(1.8)% 

NPCA 


+2.1(2.2)% 


+1.9(2.6)% 


1.455(9) 


2.427(10) 


3.532(17) 


— 


1.466(18) 


2.380(17) 


3.488(34) 


- 


+0.8(1.4)% 


-1.9(1.3)% 


-1.2(1.1)% 


- 



TABLE XL Mass splittings at ml = 2.5a 



6/g [nip - m^)/a (mjy - mA)/y/a 

WA 
5.7 1.22(3) 

6.0 1.79(2) 

FWA 
5.7 1.19(7) 0.28(8) 

6.0 1.72(6) 0.37(7) 

TFCA 
5.7 2.54(7) 0.50(10) 

6.0 2.68(6) 0.51(8) 

OFCA 
5.7 2.94(8) 0.53(10) 

6.0 2.84(6) 0.54(8) 

NPCA 
5.7 3.05(6) 0.48(5) 

6.0 3.00(2) 
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FIG. 1. Zero mode pole position expressed as a bare quark mass vs instanton size for artificial lattice instantons. Note that 
all the curves eventually approach negative infinity for small instanton sizes. 
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FIG. 2. Probability distribution of leading eigenvalue for various fermion actions on 100 12^ gauge configurations at quenched 

6/g^ = 5.85. 
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FIG. 3. Noise to signal ratio for the pion correlator at time i = 10 us m,r/n^p for the NPCA and TFCA, on a set of 80 8'^ x 24 
lattices at /3 = 5.7. 
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FIG. 4. Hadron masses in units of the string tension for various actions vs lattice spacing at fixed m^ — 2.5a. Masses are 
extrapolated to a common continuum value using a function linear in a for WA and FWA and a function finear in a^ for the 
other actions. 
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FIG. 5. Hadron mass ratios for various actions vs lattice spacing at fixed m^ = 2.5a. The extrapolation is the same as in 



Fig. 
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